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Abstract 

The Cryogenic Dark Matter Search recently announced the observation of two signal events with 
a 77% confidence level. Although statistically inconclusive, it is nevertheless suggestive. In this 
work we present a model-independent analysis on the implication of direct and indirect searches 
for the dark matter using effective operator approach. Assuming that the interactions between 
(scalar, fermion or vector) dark matter and the standard model are mediated by unknown new 
physics at the scale A, we examine various dimension-6 tree- level induced operators and constrain 
them using the current experimental data, including the WMAP data of the relic abundance, 
CDMS II direct detection of the spin-independent scattering, and indirect detection data (Fermi 
LAT cosmic gamma-ray). Finally, the LHC search is also explored. 
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I. INTRODUCTION 



The observational data that have accumulated for decades point to the existence of a 
significant amount of dark matter (DM), and the nature of this substance has become one 
of the key problems at the interface between particle physics, astrophysics and cosmology. 
Since the Standard Model (SM) of particle physics does not possess any candidate for such 
dark matter, this problem constitutes a major indication for new physics beyond the SM. 
Our current knowledge of the proprieties of DM are inferred solely from astrophysical and 
cosmological observations, which provide no information about its most basic characteristics 



such as the spin and mass. Recently. 



for example, INTEGRAL fl, ATIC 



many anoma. 



ies in various cosmic-ray observations, 
sl and Fermi LAT j4|, were found and 



2|, PAMELA 

attracted many investigators' attention. Due to our poor understanding of astrophysical 
backgrounds, it is not affirmative to say those excesses are indeed caused by dark matter. 
On the other hand, the direct detection is considered more robust. The Cryogenic Dark 
Matter Search (CDMS II) recently announced the observation of two signal events with 
a 77% confidence level ^ and additional two events just outside the signal region. While 
statistically such an observation is not significant, it is nevertheless suggestive and stimulates 
a lot of studies of its implication {g-IS]. 

^htest supersym- 



There are many proposals for the dark matter in the literature: the lig 
metric particle in the supersymmetric models with conserved -R-parity |16l . Q, the 



ightest 



KK excitation in the universal extra dimension models with conserved KK-parity 
the hghtest T-odd particle in the little Higgs models with conserved T-parity j2l|-|23|. and 



dark scalar models 



24|-l35l|. etc. Beside the DM candidate, all the new physics models 



mentioned above also provide many novel predictions that can (and will) be probed at the 
Large Hadron Collider (LHC); observing any such new effect will shed light on the theory 
underlying the SM. But it is also possible that all new physics effects will be hidden from 
direct observation at the LHC [36]. This is the case, for example, when the DM belongs to 
a "hidden sector" that interacts very weakly with the SM particles (e.g. if the mediating 
particles are very heavy), see Fig. [T] for illustration. Communication between the hidden 
sector and the SM can then be described by an effective field theory (EFT). In this paper, we 
investigate the dimension six, tree-level generated and SU{3)c x SU{2)l x U{1)y invariant 
effective operators, which describe the interactions of DM and SM particles. Similar studies 
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FIG. 1: Pictorial illustration of our EFT model, in which the SM particles interact with the 
(unknown) DM particles through the new physics which appears at the scale A. 



in the effective Lagrangian approach have been carried out in the hterature 37H46| . which 
focus on different operators. 

We add one new field D (the DM candidate) and define the effective Lagrangian as 

f(5) ^(6) 

where O-"^ is a dimension-n operator which consists of a gauge-invariant combination of the 



SM and D fields 



47 



48| : A then denotes the new physics scale (we have assumed that the 
new physics decouples in the limit A — )■ oo). Since the loop-induced operator suffers from 
l/(167r^) loop-integration factor, we focus our attention on the tree-level induced operators 
throughout this work. 

The DM candidate field D is assumed to carry an odd quantum number under a discrete 
Z2 symmetry^ under which all SM fields are even. Such a Z2 symmetry could be an unbroken 
remnant of some underlying U{1) gauge symmetry 5^. When D is a scalar or fermion, C^jf 
is then required to be invariant under SU{?>)c x SU{2)l x U{\)y x Z2. When D is a vector, 
it will be assumed to be the gauge boson of an additional U{l)x symmetry and to carry 



^ This eliminates possible contributions to observables that are strongly constrained by existing data. For 
example, if D is a vector field, the Z2 symmetry forbids contributions to 4-Fermi current-current interac- 
tions whose coefficients are tightly bound [4^. 
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no SM quantum numbers; all SM particles are assumed to be U{l)x singlets and Ceff will 
be invariant under SU{3)c x SU{2)l x f/(l)y x U{l)x- The main object of this work is to 
examine the constraints imposed by current DM searches on the coefficients We also 
explore the collider signatures at the LHC and probe these SM-DM interactions in future 
DM search experiments. 

The rest of this paper is organized as follows. In Sec. |TT] we briefly review some general 
aspects of dark matter. In Sec. Illlt HVlandrVlwe present our studies of fermionic dark matter, 
scalar dark matter and vector dark matter, respectively. Finally, in Sec. IVIt we review the 
results of our study and present our conclusions. 



II. GENERAL ASPECTS OF COLD DARK MATTER 

In order to survive until the present epoch, any dark matter particles must either be stable 
or have a lifetime longer than the present age of the universe. Furthermore, if the dark matter 



particles have electromagnetic or strong interactions, they would bind to nua 
anomalous heavy isotopes. Such isotopes have been sought but not found [51 



eons and form 



53|. Thus, the 



dark matter particles can, at best, participate in weak (and gravitational) interaction, or, 
at worst, only participate in gravitational interaction. One obvious possibility satisfying the 
foregoing constraints is that dark matter consists of neutrinos. However, the present data 
on neutrino masses show that, although neutrinos might barely account for the inferred DM 
mass density, they cannot generate the observed structure: simulations of galaxy and cluster 
formation require the cold dark matter. That is, the weakly interacting massive particles 
(WIMPs) are favored as DM candidates. 



The contribution to the energy density from a DM, Qdm, of mass rrip) is 



54| 



^ , 2 1.04 X 10^ xf 1 



Mpi y/g^ {a + Sb/xp)' 

where h = 0.73±0.04 is the scaled Hubble parameter, xp = niD/Tp with Tp being the freeze- 
out temperature, Mpi = G^^'^ = 1.22 x 10^^ GeV, and counts the number of relativistic 
degrees of freedom at the freeze-out temperature Tp. As the freeze-out temperature is 
below the DM mass roughly by a factor of 20-25, the dark matter freeze-out mechanism 
is independent of the uncertain early thermal history of the universe and possible new 
interactions at high energy scales. The parameters a and h should be derived from {crVrei), 
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the thermally averaged annihilation cross section of DM (a) times the relative velocity (vrei), 
cf. Eq. flA7|) . The relic abundance of DM in the universe has been precisely measured: 
combining the results from the WMAP Collaboration with those from the Sloan Digital Sky 



Survey gives 



55] 



flDuh^ = O.int^of^l (2alevel). (3) 

As to be shown below, the accuracy of about 10% in the expected DM abundance imposes an 
extremely strong bound on the effective operators, and leads to a strong correlation between 
the dark matter mass and the new physics scale A, e.g., the heavier the dark matter the 
larger the scale A. We focus on heavy dark matter, say mo ^ 100 GeV, all through this 
paper and will present the study of the light dark matter elsewhere. 

Besides cosmological bounds on dark matter abundance, there are other constraints on 
the properties of dark matter coming from direct detection searches of WIMPs in the halo 
of the Milky Way. The idea is that if WIMPs constitute the dominant component of the 
halo of our galaxy, it is expected that some may cross the Earth at a reasonable rate and 
can be detected by measuring the energy deposited in a low background detector through 
the scattering of a dark matter particle with nuclei of the detector. Several experiments 
have obtained upper bounds on scattering cross sections as a function of WIMP mass. The 
elastic scattering of a WIMP with nuclei can be separated into spin-independent (SI) and 
spin-dependent (SD) contributions. SI scattering can take place coherently with all of the 
nucleons in a nucleus, leading to a cross section proportional to the square of the nuclear 
mass. As a result, the current constraints on spin-independent scattering are considerably 
stronger than that on the spin-dependent component. As to be shown, the direct search of 
dark matter has a significant impact on the effective operators of both fermionic and scalar 
dark matters. 

Additional constraints on relic DM can be derived from the expectation that DM will 
collect and become gravitationally bound to the center of the galaxy, the center of the Sun 
and the center of the Earth. If this happens, then a variety of indirect dark matter detection 
opportunities arise. Detecting the annihilation products of dark matter particles in the form 
of gamma-rays, antimatters and neutrinos are collectively known as indirect search. Among 
these cosmic-ray observations, gamma-ray is thought to be more robust. Therefore, we 
focus our attention on the detection of the anomalous cosmic gamma-rays and examine the 
impact of current measurements and future projected sensitivities on the effective operators. 



Studies of antimatter and neutrino searches are also interesting but are beyond the scope of 
this project. 

If the DM particles are light enough, they will be pair-produced at Large Hadron Collider 
(LHC) and subsequently escape from the detector, resulting in a spectacular collider signa- 
ture characterized by a large missing transverse momentum. These events can be tagged by 
requiring that they also contain specific SM particles, such as hard photons, charged leptons, 
or jets. Extracting the DM signal from this type of events reduces essentially to a counting 
experiment since a resonance in the invariant mass distribution of the DM particle pair can- 
not be extracted. Thus, one needs to have better control of the SM background to affirm the 
existence of DM particles. Here we focus our attention only to the statistical uncertainty 
of the background events, but needless to say, other uncertainties, such as the systematic 
error and parton luminosity uncertainty, etc., must be included in order to provide realistic 
predictions. 

Each experiment mentioned above will have different sensitivities to each individual dark 
operator (i.e. the effective operators) and thus can probe the parameter space of the dark 
operators independently ^. Based on the time-line of the experimental programs, it can be 
achieved in the following steps: 

• Current relic abundance measurements severely constrain the dark operators and in- 
duce a non-trivial relation between the dark matter mass and the new physics scale. 

• Direct and indirect searches of dark matter, mostly the SI scattering experiment, 
further constrain the parameter space of the dark operators. 

• Dark matter pair can be copiously produced at the incoming LHC: discovering events 
with large missing energy can determine the possible mass of the DM. 

• Future direct and indirect detections of dark matter can probe more parameter space 
lying beyond the reach of LHC. 

• Consistently checking all experimental measurements may shed light on the underlying 
theory. 



Systematic studies of the neutralino dark matter along this direction have been carried out in Refs. 



56-5. 
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The main goal of the remainder of this paper is to show the close connection between dark 
matter searches and the experiments foreseen at the LHC. We will address this question 
individually for the fermion, scalar and vector dark matter below. 



III. FERMIONIC DARK MATTER x 

We first analyze the case that the DM candidate is a fermion x, which is assumed to be 
odd under a Z2, and is a SM gauge singlet; we also assume that x has no chiral interactions. 
In this case the (tree-level generated) effective interactions with the SM are of two types 

• four-fermion: 

Oex = \ {ei^e) (X7.X) , = [h] m , O,^ = (qx) (xq) , (4) 

• Vectors, fermions and scalars ^ 

O^^ = % [xipX) + h.c. . (5) 

Here, q {£) denotes the left-handed doublets while u{d, e) denotes the right-handed singlet. 
In the unitary gauge the Higgs doublet is given by 

-^(: 

where h is the SM Higgs boson and v = 246 GeV is the vacuum expectation value. All 
through this study we will choose rrih = 120 GeV. With the help of the Fierz identity, the 
scalar-scalar-current in eq (jl]) can be written in the vector- vector-current form. In this work, 
we will focus on the four-fermion effective Lagrangian for the xxff interaction in the form 
of 

4/ = (xi.Prx) {lYiglPL + 9iPR)f) , (6) 

where f = u,d,e and 

9l = -^^qx^ 9r = 2^ux, 9l = -2^qx^ 9r = 2^dx, 

9l = -\(^ix^ 9r = \oiex, 9l = -^a^x> 9r = 0. (7) 



The dimension-5 operator xx4'^4' has been studied in Refs. 59, 6^] 
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/ '\ /' '\ /'■'* ^\/^^ ^\ 

FIG. 2: Feynman diagrams for annihilations. Blobs denote the effective vertices induced by 
the dim-6 operators. 

TABLE I: Sensitivities of different experiments to the dark operators where the round brackets 
denote the Feynman diagrams in Fig. [2] while the square brackets denotes the dark operators. The 
last three columns denote the LHC collider signatures of those dark operators, 7+ (qq — ^ 7XX)) 
j+ {qq — ^ jxx) and h+ (qq — ^ ^XX); see Fig. [8] for details. 





X-nucleon 


cosmic 7-ray* 


7+ 


j+ 




(a) ] 
(d) [a^^] 


(a) [a^J^] 

(b) [a^y^] 


(a) [a^f ] 
(d) [a^y] 


(a) [ajf ] 

(b) [a^^] 


(a) [a^f ] 

(b) [a^^xl 


(d) [a^^] 

(e) [a^x] 



•k: also sensitive to the sign of the coefficients. 

Note that the right-handed projector Pr in eq originates from the fact that the fermionic 
DM particle considered in this work is a SM gauge singlet. We assume 

CXqX ~ '^/x' ^ux (-^dx ^^X ~ '^/x' 

and these operators are also diagonal in the flavor space for simplicity. The fermionic DM 
can also annihilate into a vector-scalar boson pair via the following vertex: 

4'i = ^«<^x ixYx) Z,h. (8) 
The operator Ocf,^ also gives rise to the vertex Zxx as 

41 - ^^A-^'X. (9) 

which contributes to the annihilation processes of XX Z ^ ff, Zh or W~^W~. For 
simplicity we consider the operators Cj^, Of^ and O^^ and use the coefficients, a^^^ and 
a^^, to denote the operators hereafter. 

The fermionic DM annihilates into the SM fermions, vector bosons, or into a pair of 
Higgs boson and Z-boson. The Feynman diagrams of the annihilations are shown in Fig. [2l 
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where the diagram (a) is related to and the diagrams (b-e) are related to a,^^. These 

five diagrams have different effects on various detection experiments; see Table HI For the 
annihilation of a heavy x, the contributions from diagrams (b), (c) and (e), being s-channel 
processes, are highly suppressed by l/(4m^) from the s-channel propagator. Therefore, 
only diagrams (a) and (d) need to be considered for a heavy x annihilation. It is worthy 
mentioning that if x is very light, say ~ mz/^, then the diagram (b) receives a large 
enhancement from the threshold and becomes dominant. Detailed study of such light dark 
matter will be presented in a future work. For elastic scattering from a nucleus, only 
diagrams (a) and (b) contribute, whereas for cosmic gamma-ray detection, only diagrams 
(a) and (d) contribute sizably. Moreover, diagrams (a, b) can be probed at the LHC using 
the signature of mono-photon plus missing energy while diagrams (d, e) can be probed using 
the signature of single Higgs scalar plus missing energy, and diagram (c) can be probed using 
the signal of two forward jets with missing energy. All these issues will be explored in details 
in the rest of this section. 

A. Relic abundance 

One can calculate the cross sections of XX annihilation and obtain the leading terms (a 
and h) in the non-relativistic expansion. The xp can be determined using Eq. ( 1A8I) . After 
substituting xp and (a, h) into Eq. ([2]), one can evaluate the relic abundance of the DM. 
Detailed calculations are given in Appendix [B] 

The allowed parameter set {m-^, A) with respect to WMAP data is shown in Fig. |3l For a 
heavy say > 300 GeV, all the SM particles can be treated as massless. One can then 
obtain the leading terms a and h as follows: 

2 

TfL 

a = ^ (0-24a" + 0-21"/x + O-OSaJJ , (10) 

9 

Tfl 

b = ^ (0-ll«?x' + 0-10«/x + • (11) 

We note that xp hardly varies, xp ~ 25 ^, in the entire allowed parameter space, but the 
freeze-out temperature Tp varies from 20 — 80 GeV for ~ 500 — 2000 GeV. Choosing 

^ In the equation that determines Tp, x appears in a Boltzmann factor . Taking the logarithm, one can 
show Xi? w 25 for a wide range of values of the annihilation cross section. 
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FIG. 3: Allowed parameter set (m-^, A) for a universal coupling constant a when DM is a fermion. 
The upper (lower) boundary of each band corresponds to the upper (lower) limit of O^M^^- 

xp = 25, we obtain the following relation between and A, 

2.22 X 10-2 _ / \2 /TeV 

(nDMh^) (0.252a^^2 + 0.221a^^2 ^ o.OSlaJ^) ~ UOO GevJ 

We have verified that this relation agrees very well with the numerical results shown in Fig. 

El 

For illustrations, we adopt the value of a to be 1 in the following discussions, but it 
is straightforward to include its effects by rescaling m. When = 500 GeV, the second 
lightest new particle is expected to have a mass of about 2 TeV if only aj^ i^fx) non-zero, 
but it should appear around 1 TeV if a,^^ is present only, see respectively the red (green) and 
blue bands. The bands of aj^ (red) and af^ (green) overlap while the aj^ band is slightly 
higher because the coefficient of a^^ is slightly larger than the one of a^^, see Eq. f|T2l) . This 
small difference comes from three neutrino annihilation channels which only contribute to 
aj^ but not to af^. Furthermore, the coefficients of a^^ and af^ are both larger than the 
coefficients of a^,^ by an order of magnitude, i.e., the annihilation into fermions is dominant 
in the annihilation cross sections. Eq. (jj]) implies that, for a given m, the scale A has to be 
small in order to compensate the smaller coefficient of a^^, leading to a much lower band. 
Finally, when these three couplings all are present, the allowed new physics scale A becomes 
larger, e.g. A ~ 2.6 - 5.5 TeV for ~ 500 - 2000 GeV. 

It is worthy mentioning that for each operator the region below the band is also allowed 
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even though less rehc abundance is produced and additional annihilation channels are needed 
to explain the current relic abundance. 



B. Direct detection of x 

The elastic scattering of a fermion dark matter x from a nucleus is induced by effective 
four fermion operators: 

^ f ^.R ^ , ^ 9 V 



+ ^2 + + "^x^^^— ) X7.75XgV75g, (13) 

with Cy = Tl — 2Qg sin^ 9]v and C\ = T|. Here, T3 and Q are the weak isospin and electric 
charge of the quark q, respectively, and 6w is the Weinberg angle. The first term results in a 
spin-independent scattering from a nucleus while the second term leads to a spin-dependent 
scattering from a nucleus. It is convenient to consider the cross section with a single nucleon 
for comparing with experiments. 

The DM-nucleus cross section in this case is ^ 

SI _ ^I'^nK ,x 

where is the effective DM-nucleus coupling. Through simple algebra (see Appendix IC II 
for details), we obtain the DM-proton and DM-neutron cross sections as given below: 

(TeV \ ^ 
-^j «-«J^ + 0.013a<^^)', (15) 

a^i ^ (6.98 X 10-^ pb) (^^j « - a^^ - 0.162a^J , (16) 

where ^ nip ~ m„ are used to derive the approximate results. We note that cr^p ~ a^^ 
when either a^^ or is non-zero, whereas ^ cr^^ when only presents. 
The spin- dependent DM-nuclei elastic scattering cross section can be expressed as 



a 



sn^ /'""f" [Ay(J+l)]. (17) 



^ Note that Eq. ((Ti| does not correspond to the total DM-nucleus cross section, but to that at zero mo- 
mentum transfer, usually named as "standard cross section" . 
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where J is the total angular momentum of the nucleus and Aat (oc 1/ J) depends on the axial 
couplings of DM to the quarks. For odd-proton nuclei the spin-dependent DM-nucleus cross 
section is mainly due to the DM-proton interactions, whereas for odd-neutron nuclei it is 
dominated by DM-neutron scattering. For even-even nuclei the spin-dependent cross-section 
is highly suppressed. For the proton/neutron as the target, Eq. fll7p is transformed into the 
cross section from DM-proton/neutron interactions with the proton/neutron spin, which are 
given as 

SD 

XV 



(af^ + 05-^-4.530^;,)', (18) 
^ (4.183 X 10-6 pb) {af^ + a^^ + 3.53a,,)' . (19) 

It indicates that for a given A the contribution of O^^ dominates over the other two operators. 
Furthermore, we note that a^E' ~ cr^f* 

' xy X"' 

In Figs. Ht^a) and|4](b), we compare the spin- independent elastic scattering cross section 
of the fermionic DM x with the current and projected sensitivities of direct detection exper- 



iments; the blue-solid curve denotes the XENONIO exclusion limit 
shows the CDMS II 2009 exclusion limit 



6l| while the red curve 



62| . In the near future, the CDMS II is projected 



to increase the sensitivities of spin-independent scattering roughl y by a factor of 4; see the 
black-dashed curve. The projected sensitivities of SuperCDMS [63|, whose first phase is 



proposed to start operating in 2011, are also plotted; the blue-dashed line represents phase- 
A while the red-dashed line phase-C. The red, green and blue bands represent the cross 
sections for a^^ = 1, af^ = 1 and = 1, respectively, with the assumption of only one 
parameter being non-zero at a time. We note that the bands of a^^ and aj?^ overlap for the 
whole mass range we are interested, whereas the band of a^^ is much higher than those two. 
It is consistent with the relic abundance constraints shown in Fig. |3l The exclusion limit 
of (denoted by the symbol ">") and the potential reach of future experiments (denoted 
by "<") are summarized in Table HTl If the recent CDMS II observation is a hint that di- 
rect detection is "around the corner" , the the dark matter mass should hide in between the 
CDMS II 2009 exclusion limit and CDMS II projected sensitivity. Nevertheless, the entire 
range of considered here can be covered by the Phase-C of SuperCDMS. 

Considering now the case when all three operators contribute, for which the scattering 
cross sections are shown as the broad gray bands. Although the relic abundance imposes very 
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[GeV] [GeV] 

FIG. 4: Prediction of the fermionic DM- nucleoli cross sections with respect to m-^ and A for 
aUowed parameter set given in Fig. [2 (a, b) for the spin-independent DM search; (c, d) for the 
spin-dependent DM search. The red (green, blue) band denotes the cross section for = 1 
^<t>x ~ -*■)' respectively, when one operator is considered at a time. The gray shaded 
region denotes the cross sections when all these three operators contribute, which corresponds to 
the black band in Fig. [3l see the text for details. In the upper panel, the blue-solid (red-solid) 
line labels the upper limit (90% confidence level) on the SI DM-nucleon scattering cross section 
from the current XENONIO (CDMS II 2009) direct search, respectively. The black-dashed line 
denotes the near term expected sensitivity from the CDMS II experiment, whereas the blue-dotted 
(red-dotted) line represents the longer term projection for the Super-CDMS Phase-A (Phase-C), 
respectively. 
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TABLE II: Exclusion limits of (GeV) from current SI direct search (denoted by the symbol 
">") and potential reach of future experiments (denoted by the symbol "<") when one parameter 
is considered at a time. The color of each operator refers to Fig. HI 

XENONIO CDMS II 2009 CDMS II proj SuperCDMS-A SuperCDMS-C 
Q^^ (red) > 750 > 1050 < 2000 < 2000 < 2000 

af^ (green) > 800 > 1130 < 2000 < 2000 < 2000 

a0^(blue) > 600 > 870 < 1700 < 2000 < 2000 



tight constraints on m^, A and the coefficients, see Eq. (IT^ . the signs of the coefficients are 
not determined as the rehc abundance is only proportional to the square of the coefficients. 
On the contrary, the scattering cross sections are sensitive to the signs of the coefficients; see 
Eqs. fllStiigp . The upper limit of the x-proton SI scattering and the neutron SI scattering 
corresponds to a^^ = —df^ = «0x ~ ~^fx ~ ^fx ~ ~ respectively. The 

lower limits of both scattering corresponds to aj^ = aj-^, i.e., solely determined by the 
contribution of a^^. The large difference between two lower limits, for which x-neutron 
scattering is much larger than x-proton scattering, is attributed to the the coefficients of 
a^^ in Eqs. (fTS!) and (fT6|) . 

For completeness, in Figs. |l](c) and |l](d) we also present spin-dependent cross section 
which is obviously well below the present and even future experimental reach. Picking 
out such a signal amongst the large range of noise sources is a considerable challenge. See 
Ref. iJ] for a brief comment on the impact of spin-dependent measurement on the PAMELA 
antiproton observation. 



C. Indirect gamma-ray detection of x 

Photons from dark matter annihilation in the center of the galaxy also provide an indirect 
signal. Detection of cosmic gamma-rays is carried out in both space-based and ground-based 
telescopes; the former directly observe the cosmic gamma-rays and can cover the low energy 
regime (GeV), whereas the latter indirectly observe the gamma-rays through the detection 
of secondary particles and the Cerenkov light originating from their passage through the 
Earth's atmosphere and will cover the high energy regime (100 GeV to TeV). Thus, both 
types of telescopes can be complementary to each other. The sensitivity of present and 
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Log[ E [GeV] ] Log[ E [GeV] ] Log[ E [GeV] ] 



FIG. 5: (a) Predicted gamma-ray spectra for the annihilation of fermionic DM x (sohd hnes) with 
the NFW density profile. The Fermi LAT observation (Galactic background) is also presented by 
the red-box (grey-solid line); (b) Comparison between the DM signal plus background and the 
Fermi LAT observation. Note that a common boost factor of 20 is applied to all the DM signals; 
(c) The comparison between the DM signal with the difference of Fermi LAT observation and 
background. The different boost factors are adopted for different DM masses so that the DM 
signal will not exceed the data. 



future detectors in gamma-ray astrophysics has been studied in Ref. 64 1. 

As the DM x cannot annihilate into photon pairs directly, we can only detect continuum 
photon signals. Using Eq. ( ID2|) we obtain the differential flux of the gamma-rays observed 
on Earth from DM annihilation as follows: 

'^^ i (1.4 X lO-^^s-^^cm-^GeV-i) J(An) Afi 



X (^^) (^) + "fx + 0.042e-"-a5J , (20) 

where x = E^/ and J (AQ) AQ counts the dependence on the DM halo profile (see 
Appendix [D] for details). For simplicity, we assume a standard NFW density profile 65l. |66| 
for the DM in our galaxy, i.e. J (AQ) AQ ~ 1 for a AQ = 10^^ sr region around the direction 
of the galactic center. The predicted differential gamma-ray fluxes from DM annihilation 



are plotted in Fig. [5]^a) with the existing Fermi LAT observations (red box) 67|, |68| and the 



galactic background (gray-solid curve). See Ref. [69| for details of the galactic background 
and its uncertainties. For illustration, we choose as 30 (50, 100, 300, 1000) GeV with 
the corresponding scale A which is derived for each mass from the relic abundance, cf. Fig. 
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[3l We only present the distribution of a^^ = 1 in the figure as all three operators give rise 
to almost the same distributions. It is because the shape of gamma-ray spectrum is almost 
independent of the specific annihilation process and the normalization is also fixed by the 
relic abundance. Clearly, all fiux distributions of the gamma-rays are consistent with the 
current Fermi LAT observations since the dark matter signals are far below the observation 
and background. 

However, the distribution of dark matter might be clumpy. Such small-scale structure 
would enhance DM signal by a boost factor (B) defined in Eq. (]D3I) . Due to large uncer- 
tainties of Galactic background estimations, even there is small discrepancy between signal 
and background, as shown in Fig. O one should not treat such small deviation too seriously. 
Instead, one can extract upper bounds on the allowed boost factor for dark matter. We 
vary the boost factor such that the distributions of cosmic gamma-rays for the given 
plus background satisfy the Fermi LAT measurement. When the coefficients of the dark 
operators are set to be 1, we found that a boost factor S = 20 is allowed by the Fermi LAT 
data for the 30 GeV dark matter; see Figs. Et^b) and (c). However, it is in contradiction to 
the CDMS II observation when one interprets the two events as signal. For a heavier dark 
matter (m^ > 30 GeV), the allowed boost factor can be larger. With a large boost factor, 
such a heavy dark matter would be detected by Fermi LAT with good sensitivity in higher 
energy region since the dark matter annihilation produces much harder energy spectrum of 
the cosmic gamma-ray; see Fig. |5t^c). 

The integrated photon fiux is also interesting. The integrated photon fiux above some 
photon energy threshold Eth is given by Eq. ( 1D4I) . We plot the integrated photon fiuxes 
in Fig. [6] as a function of for two representative Efh- 1 GeV, accessible to space-based 



detectors, and 50 GeV, characteristic of ground-based telescopes, 
two promising experiments, Fermi LAT TOj and MAGIC II 71 



jstimated sensitivities for 



72| . are also shown. It is 



clear to see from the Fig. [6] that the Fermi LAT and MAGIC II could probe a heavy x if 
the background is well understood. 



D. Collider search for x 

Since xp ^ 25 a.s pointed out above, one can substitute it into the equation of the relic 
abundance Eq. ([2]) and obtain the proper thermally averaged annihilation cross section. 
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(GeV) 

FIG. 6: Integrated photon flux as a function of for energy threshold of 1 GeV (upper bands) 
and 50 GeV (lower bands). The notation of the color band is the same as Fig. H]). The plot assumes 
J{^, Af2)AQ = 1; all fluxes scale linearly with this parameter, see the Appendix ID] for details. 

which gives rise to the correct relic abundance, as {3] 

{(^annVrel) ~ 1 pb. 

Inverting the annihilation process, we can produce the DM pair at the LHC though the 
initial state only consists of light quarks. The XX P^^^ can be produced via the processes 

qq XX, qq-^ Z ^ XX, (21) 

Shown in Fig. [7] is the cross section of the xX P^ir production (aprod) as a function of 
m^. The upper narrow band denotes the cross section for aj^ = 1 (red), or af^ = 1 
(green), or aj^ = af^ = a^y^ = 1 (black), whereas the blue narrow band below denotes 
the cross section for a^^ = 1. The upper three narrow bands almost degenerate due to 
the relic abundance constraints. We note that <Jprod ~ 1 pb for ~ 100 GeV, but <7prod 
drops steeply with increasing m^. The large suppression is due to two facts. The first is 
the decrease of coupling, since the A gets larger. The second reason is due to the parton 
distribution functions (PDFs) which drop rapidly in the large Bjorken-x region. Typically, 
(x) ~ 2my,/y/s with y/s = 14TeV for LHC and y/s = 1.96 TeV for the Tevatron. Such a 
large suppression restricts the DM search at the LHC to small m^. The a<^^ operator suffers 
from a much severer suppression as it contributes only via the Z-boson-exchange process, i.e. 
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[GeV] 

FIG. 7: Cross section of the XX P^^^ production at the LHC. The lower narrow band (blue) denotes 
the cross section for = 1 while the upper narrow bands (black, green and red) denote the cross 
sections for a^^ = 1, aj^ = 1, or a^^ = aj^ = a^^ = 1. The shaded regions denote the preferred 
DM mass region by the CDMS II positive signal; see Table HIl 

the second term in Eq. fl2T]) . For a heavy DM pair, the off-shell Z boson in the propagator 
receives a large 1/s suppression which makes the cross section drop much faster than other 
bands. The shaded blue (green, red, gray) region denotes the preferred dark matter region 

for a^^ ('^/x' '^/x' '^/x ~ ^fx ~ ^'i'x ~ ^) CDMS II positive signal. 

In order to detect the DM signal, additional SM particles are needed. A hard photon 
is a good probe which can be produced in association with the DM pair via the following 

processes 

QQ qq-^ Z'y XX7, (22) 

as shown in Fig. [8](a) and (b), giving rise to a collider signature of mono-photon plus large 
missing transverse momentum (7+ ^t)- Similarly, one can also search for the DM x using 
the signature of mono-jet plus missing transverse momentum (j+ ^t) which involves the 
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jtC^x y^v#i^ 

(a) (6) (c) 




FIG. 8: Schematic diagrams of the searching channels of DM x the cohiders. 
processes as follows: 

QQ XX9, Zg ^ XX9, 

qg xxq, qg Zq^ xxq- (23) 

Furthermore, the XX P^ir can also be produced associated with a Higgs boson and through 
the W-boson fusion process, i.e. 



qq 



Z-^XXK qq-^Zh^xxK qq' ^ W+*W *qq' ^ Z*qq' ^ xxqq' , (24) 



as shown in Fig. [Ht^c), (d) and (e). Unfortunately, the cross sections of such processes Eq. 
(El!) too small to be detected. We then focus our attention on the 7+ and j+ 
signatures below. 

Search for new physics using the collider signatures of 7+ I^t and j+ I^t has been carried 
out at the Tevatron Run-I iJl = 1.8 TeV) and Run II (\/s = 1.96 TeV ), which set upper 

nn 

bounds on the production cross sections of above processes [7J, 175[. For the parameter space 
given in Fig. |3]the production cross section is below the Tevatron direct search bound and is 
not shown here. We plot the relevant cross sections at the LHC in Fig. |9l (a) 7+ and (b) 
j+ ]^T- In Fig. IHl^b) we have summed over the processes shown in Eq. ( l23l) . Similar to Fig. [TJ 
the lower blue narrow band denotes the cross section for a^-^ = 1 while the upper narrow 
bands denote the cross sections for a^^ = 1 (red), af^ = 1 (green) and a^^ = af^ = a^j,^ = 1 
(black), respectively. In order to avoid the collinear singularities arising from the light quark 
propagator, we impose a kinematical cut on the transverse momentum (px) of the final state 
photon or jet as pp'''^ > 5 GeV. The shaded regions denote the preferred DM mass region 
by the CDMS II positive signal (see Table HTll . 

Besides the PDF suppression, the production of the DM pair associated with one single 
photon (or with one single jet) is suppressed by two additional factors: (i) additional coupling 
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FIG. 9: Production cross sections (in the unit of fb) of the x P^'ii' together with a SM particle at 
the LHC: (a) mono-photon plus ^t-, (b) mono-jet plus ^t-, where a kinematics cut > 5 GeV (or 
p?P > 5 GeV) is imposed in order to avoid the collinear singularities arising from the light quark 
propagators. The lower band (blue) denotes the cross section for a^^. = 1 while the upper bands 
(black, green and red) denote the cross sections for a^^ = 1, or aj?^ = 1, or a^^ = aj^ = a^-^ = 1. 
The shaded regions denote the preferred DM mass region by the CDMS II positive signal; see 
Table IE 

{aEM or as) and (ii) large volume of three-body phase space. We note that the cross sections 
of the mono-photon production are much smaller than the ones of the mono-jet production. 
It is because the former only comes from the quark-quark initial state while the latter comes 
from both the quark-quark initial state and the quark-gluon initial state. Although having 
a smaller cross section, the 7-I- I^t process has a relatively clean signature, i.e. hard photon, 
which could be detected at the LHC. On the other hand, the mono-jet process suffers from 
huge QCD backgrounds which makes its detection more challenging. Furthermore, the PDFs 
receive large uncertainties in the large DM mass region, i.e., the large x region. Such large 
uncertainties will further make the detection of the signal event intricate. 

The blue narrow band, i.e. the cross section for a^^^ = 1, is much lower than other three 
narrow bands because a^^ only contributes via the Feynman diagram (b) in Fig. [HI Bearing 
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in mind that Higgs-associated production is also rare, we then conclude that it is hard to 
observe the sole a^y^ contribution at the LHC, i.e. one can probe a^-^ only via the cosmology 
observations. For the other operators, the positive CDMS II signal prefers the DM mass 
~ GeV, where the mono-photon and mono-jet production cross sections are of 
order (9(10~^) fb and 0{1) fb, respectively. It is plausible to cross-check the cosmological 



and collider measurements when the good knowledge of the SM 



Dackgrounds and detector 



sensitivity at the LHC is achieved. For example, see Refs. |40|, |42|, |45| for a detailed study of 
bounds from existing Tevatron searches for monojets as well as expected LHC reaches for a 
discovery. 



IV. SCALAR DM Lp 

Let's now consider the scalar type dark matter. After adding a Z2-odd real scalar singlet 
y9 into the SM, we obtain the effective operators as follows: 



scalar only 



scalar-vector 



scalar fermion 



O^i = I (0+0) ' ivv) , = d, (0+0) (^^) . (25) 



0^3 = (ipip) {D,(l>+D^(f) (26) 



= (v^V^) (^e0) , Oucj, = {(f(f) [qu4>^ 

Od^ = (v^v^) {qdct>) . (27) 



Besides the operators listed above, there exists dimension-4 operator {fLpcjycj) and dimension 



5 operator (fifff which has been studied in Ref. j60|]. There are also a few loop-induced 
operators which are suppressed by factors of 1/(167?^). In this work we only consider the tree- 
level induced operator, but it is worthy mentioning that those loop-induced operators, e.g. 
ifipB^'^ B^y, are also very interesting as they may provide a spectral line feature observable 
in indirect detection experiments if their contributions dominate in dark matter annihilation 

3- 
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These operators induce the following DM annihilation processes, 



—J- hh, ifip — > VV, and (p^p — )■ //, 

where /i(V, /) denotes the SM Higgs boson (gauge boson, fermion). The vertex of the DM 
annihilating into the Higgs scalars is 

( f\\ 1 

^)p!phh = j^9whh VVhh, (28) 

where 

g^^hh = —Oi<i,i + - (4a02 + a^a) (29) 

with s being the square of the system energy. The vertices of the DM annihilating into the 
vector bosons are summarized as follows, 

^tlvv = ("<^3 cpcpW^W-^ + ^a^3 m| ifcpZ.Z^'^ . (30) 
The vertices of the DM annihilating into the fermions are 

C'ff = -^^^ (^-^ae^CLCR + -^ttu^ULUR + -^ad^dLd^^ . (31) 
In this work we consider two possible scenarios of the (f^pff operators: 

(A) Universal a/^, i.e. = au(f, = Oidcj,. As to be shown later, such a case is ruled out by 
the current DM direct search experiments. 

(B) The coefficients are proportional to the fermion mass, i.e. aj^v = rrifa'j^, and a'j-^ is 
universal for leptons and quarks, i.e. a'^^ = a'^^ = a'^^. For simplicity, we will use a/^ 
instead of a'^^ hereafter. 

Finally, the operator C^i induces the vertex ipiph 

= '^T^^l^' (32) 

which contributes to the processes of DM annihilation into the SM Higgs bosons, vector 
bosons and fermions through the Higgs-mediated s-channel processes. In summary, there 
are four operators contributing to the scalar DM ip annihilation: O^pi, 0^2, and Of^. 

When the scalar ip is the DM, we must consider a large number of annihilation processes, 
e.g. ff/WW/ZZ/hh and ipip h ^ ff/WW/ZZ/hh, see Fig. M Sensitivities of 
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FIG. 10: Feynman diagrams for (pip annihilation. Blobs denote the effective vertices induced by 
the dim-6 operators. 



TABLE III: Sensitivities of various experiments to the dark operators where the round brackets 
denote the Feynman diagrams in Fig. [TD] while the square brackets denotes the dark operators 
given in Eqs. ([251132]) . 





X-nucleon 


cosmic 7-ray 


7+ 




VBF 


(a) [a/0] 


(a) [oLf^] 


(a) 


(a) [a/0] 


(a) [a/0] 


(c) [a^a] 


(c) [oc/ia] 


(b) [a^i] 








(d) [a0i] 


(e) [a0i,</,2,03] 




(e) [a0i,02,<^3] 









these diagrams on various experiments are summarized in Table IIIII The contributions from 
the diagrams (b,d,f) are large only around the Higgs resonance region, i.e. ~ ?^/i/2, but 
for the annihilation of heavy if pair their contributions are highly suppressed. As we are 
interested in the region of > 100 GeV, only the diagrams (a,c,e) need to be considered. 
For the elastic scattering from the nucleus, only diagrams (a) and (b) contribute, whereas 
for the cosmic gamma-ray detection, only diagrams (a,c,e) contribute. Moreover, diagrams 
(a,b) can be probed at the LHC using the signature of mono-photon plus missing energy 
while diagrams (c,d) can be probed in the vector-boson-fusion process qq — )■ VVqq — )■ ffqq. 
Diagrams (e,f) cannot be probed at the LHC. 
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FIG. 11: (a) Allowed parameter set (m,^, A) for the scalar DM ip: (a) scenario-A (family universal 
coupling); (b) scenario-B (family dependent coupling). The gray shaded region is excluded as 
rriip > A. The upper (lower) boundary of each band corresponds to the upper (lower) limit of 
^Duh^ given in Eq. ([3]). 

A. Relic abundance 



The scalar DM can annihilate into light fermion pairs, or WW and ZZ pairs if the 
vector boson channels open. The calculation of the relic abundance is similar to the one 
of fermion DM x- We present the annihilation cross sections and the leading terms in the 
non-relativistic expansion in Appendix [Bl The allowed parameter space of (m^, A) is shown 
in Fig. [TT]for both scenario-A and scenario-B. 

If the DM annihilation is purely induced by the operator O^i, the scale A decreases 
rapidly with increasing in order to respect the WMAP data; see the red band. It can be 
easily understood from dimensional counting as follows. For a heavy niip the effective vertices 
induced by C</,i give rise to the leading terms a and h as following; see Eqs. ( IB18tlB22|) . 

a, b (X — TTTT- (33) 
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Since a and b are fixed by the well-measured relic abundance, 

A « m-'/\ (34) 

which results in the inverse behavior of the red band. Moreover, is about equal to 
the cutoff scale A around 500 GeV. We thus conclude that if is the only source for 
DM annihilation, then ip should be very light, say m^p < 500 GeV. Otherwise, other dark 
operators have to be considered to explain the relic abundance. 

The bound on the operator Oftf, is not sensitive to m^; see the magenta bands. It can 
also be understood from dimensional counting. For a heavy (p, the effective vertices induced 
by Of^ give rise to the following leading terms; see Eqs. (IB13tlB16|) : 

— , (scenario — Aj, 

a, 6 oc K (35) 
\^—, (scenario -B), 

therefore, the relic abundance only depends on the scale A for both scenarios. We found 
that for all the mass region concerned, A ~ 2TeV in the scenario- A, whereas A ~ 1.1 TeV 
for scenario-B. Again in order to satisfy the condition < A, ip cannot be very heavy. For 
example, should be less than 1 TeV for scenario-B. 

Let's consider now both Ofj,2 and O^^ simultaneously, an interference between them oc- 
curs. From the annihilation cross sections given in Eqs. flB17|l and flB20|) . we derive the 
coefficients in the thermal average taking m ^ v and find 



ml 

a — 



647rA4 ^^^^i^ + 8a<^2a<^3 + ISaJg) , (36) 

9 

TTl 

25gJ'^4 (^Sajs + 8a^2a^3 + Tajg) . (37) 



The allowed parameter set (m^. A) for each operator is shown in Fig. [TTl see the green 
band (0:02) and the blue band {a^f^^). Similar to the case of fermionic DM, we find that xp 
barely varies in the entire allowed parameter space, xp — 25, but the freeze-out temperature 
Tp varies from 20 — 80 GeV for ~ 500 — 2000 GeV. Choosing xp = 25, we obtain the 
following relation between and A, 

6.98 X 10-^ _/ y/TeVV 

{ncDMh^){0.273al^ + 0.13Qa^2a4,3 + 0.213al^) VlOOGeVy \ A J ' ^ ^ 

Finally, we consider the case of all the operators contributing to the DM annihilation; see 
the black band. It is clear that for a heavy ip the DM annihilation is dominated by (9<^2 and 
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given by Eq. (138|) . Note that the sign of the coefficients a,f,2 and a^s is very important. 
For simphcity, we only consider the case of constructive interference, i.e., both a^2 and q;</,3 
being positive, in this study. But one can easily get the result of destructive interference 
by lowering the black band by ~ 25%. For a light (p, however, the DM annihilation is 
dominated by Of^ in scenario-A while by 0(f,i in scenario-B. 



B. Direct Search 



When the DM is a scalar (p, it can be detected in the spin-independent experiments via 
the scalar interaction with nucleus 

C^^m = J^^^^ll + A^^";;^^^^^' ^^^^ 
which leads to a total spin-independent <y9-nucleon cross section as follows: 

2 2 

where the effective v?-nucleon couplings, /^p and f!^, are given in Appendix IC 2[ In practice, 
we found a^j, ~ a^^. 

The SI elastic yj-nucleon scattering cross sections are plotted in Fig. [12] along with the 
current limits placed by the CDMS II 2009 (red-sohd) and XENONIO (blue-solid) collabora- 
tions, the projected CDMS II sensitivity (black-dashed), and the projected future sensitivity 
of Super CDMS, Stage- A (blue-dotted) and Stage-C (black-dotted). If the scalar DM anni- 
hilation occurs purely through the 0(f)i operator (see the red band), then the current direct 
detection from CDMS and XENON has no constraint on the DM as the elastic scattering 
cross section is far below the current bounds, almost by two order of magnitude. However, 
the operator can be probed completely in the projected future SuperCDMS at Stage-C. One 
should keep in mind that the mass of ip should be less than 500 GeV, otherwise the mass of 
if is larger than the NP scale A. However, the positive CDMS II would rule out O^i when 
the DM signal is observed in the near future. 

The magenta band shows the cross section for aj^ = 1. Consider scenario-A first. The 
cross section is so large that the entire parameter space 100 GeV < < 2 TeV is excluded 
by the CDMS and XENON data. If all the operators contribute to the DM annihilation, 
the relic abundance allowed scale A becomes large so that the SI scattering cross section is 
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FIG. 12: SI nucleoli scattering cross sections as a function of m^. The red band represents the 
cross section for q^i = 1, the magenta band denotes the cross section for aj^ = 1, and the gray 
band denotes the cross section for «/</, = a<^i = a^2 = 0(413 = 1- See Fig. [Hfor description of other 
curves. 

suppressed, see the gray band. The CDMS and XENON measurements excluded the region 
of < 1450 GeV, and the projected CDMS II sensitivity covers up to 2000 GeV. On 
the contrary, in scenario-B, the SI scattering cross section is well below the current CDMS 
and XENON data. Even the long-term SuperCDMS at Stage-C can only probe up to 
300 GeV. Similar to the case of fermion DM x, the broad gray band is due to the different 
sign assignments of aj^ and a^i. Hence, we find following interesting points: 

1. When the CDMS II observe dark matter signal in the near future, the dark matter 
in the scenario- A (family universal coupling) lie between the mass window (1450 GeV 
- 2000 GeV) when = a^i = a(/,2 = <y<i,3 = 1- The dark matter in the scenario-B 
(family dependent coupling), however, is very light, say 100 GeV < < 150 GeV 
with A ~ TeV. 

2. When no dark matter signal is observed within the projected CDMS II sensitivity, the 
scenario-A is almost excluded for < 2 TeV. However, most of the parameter space 
of scenario-B is still allowed and is hard to probe even with the long-term SuperCDMS 
at Stage-C. 
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FIG. 13: Predicted gamma-ray spectra for the scalar DM ip with the NFW density profile: (a, b, c) 
"01 = 1, (d, e, f) aj^ = 1, and (g, h, i) a^2 =1 or a^s = 1 or a^i = a^2 = a</.3 = = 1- The 
left panel shows the predicted gamma-ray spectra for the scalar DM ip. The middle panel shows 
the the summed gamma-ray spectra of signal plus background with the Fermi LAT observation. 
Note that a common boost factor is applied to all the DM signals. The right panel shows the DM 
signal with the difference of Fermi LAT observation and background. Different boost factors are 
adopted for different DM masses so that the DM signal will not exceed the data. 

C. Indirect search 

We now consider the indirect search of ip via cosmic gamma-rays. As most of the pa- 
rameter space of scenario-A is excluded, we focus on scenario-B hereafter. Similar to the 
fermionic dark matter in previous section, we plot the differential distribution of the gamma- 
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(GeV) nip (GeV) 

FIG. 14: Integrated photon flux as a function of m^p for energy threshold of 1 GeV (a) and 50 GeV 
(b) (notation of the color bands is the same as Fig. fTTI) . The plot assume J (^', Ail) Afi = 1; all 
fluxes scale linearly with this parameter. 

ray flux for various scalar DM masses and boost factors in Fig. [131 Figs, (a, b, c) shows the 
distributions of a,^! = 1; Figs, (d, e, f) display the distributions of a/^ = 1; Figs, (g, h, i) 
denote the distributions of 0,^2 = 1, or a^s = 1, or a^i = a^2 = a^s = a/,/, = 1 as their 
distributions are almost the same. It is clear that the allowed boost factors for the operator 
afij) is much smaller than that for a^i and 0^2 (or a^^ = 1, or a^i = a^2 = Ci<j>3 = «/</> = !)• 
The integrated cosmic photon-flux as a function of mo is plotted in Fig. [TH We choose 
two energy thresholds Eth = 1 GeV and Eth = 50 GeV to mimic the space-based and ground- 
based telescopes. For simplicity we assume J {AQ) AQ = 1 and all the fluxes scale linearly 
to account for other DM density proflles. Fermi LAT and MAGIC have great potential to 
observe anomalous gamma-rays. With an enhancement factor, say ~ 10, from either the DM 
density proflle or the boost factor, almost all the allowed mass region (100 GeV — 2000 GeV) 
can be explored by Fermi LAT and MAGIC. 
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FIG. 15: Schematic diagrams of the searching channels of DM x t^is coUiders. 
D. Collider search 

Now consider the direct search of at the LHC The scalar DM ip can be detected in the 
following processes 

qq 7(^(^, (41) 
qq VVq'q 44 (42) 

where the former is from the four-particle interaction; see Fig. [T5](a). while the latter is from 
the so-called vector boson fusion (VBF) process; see Fig. [TST b) and (c). 

The collider signature of the process in Eq. (14T1) consists of a hard photon plus large 
fiT due to the missing DM particles. Furthermore, can also be probed with a collider 
signature of a mono-jet plus large The corresponding processes are as follows: 

qq gpp, qg qipip. 

The signal of the VBF process is characterized by two quark jets, which typically stay in the 
forward and backward regions of the detector and are widely separated in pseudo-rapidity, 
and also by a large missing transverse momentum {^t), due to the two missing DM particles. 

In Fig. [16] we plot the cross sections of the above processes as a function of with the 
corresponding A consistent with the DM relic abundance. As only O^i and Of^ operators 
can contribute to the signatures of 7+ and j+ ^t, we plot their effects separately in 
Fig. [16] (a) and (b) where the red curve denotes 0(f,i while the magenta curve denotes Of,p. 
Although both O^^ and are involved in the VBF process, the latter contributes much 
less than the former, especially for a heavy p. Therefore, we only present the cross section 
of 0^3 in Fig. [TBT c). We note that the cross sections are generally small for all the processes. 
Note also that various kinematics cuts are needed in order to suppress the SM background. 
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FIG. 16: Production cross sections (in the unit of fb) of the if pair together with a SM particle 
at the LHC: (a) mono-photon plus fJx, (b) mono-jet plus fJx, and (c) VBF. Here, a kinematics 
cut pj,>5 GeV (or pj.>5 GeV) is imposed in order to avoid the collinear singularities arising 
from the light quark propagators. In (a) and (b), the lower band (red) denotes the cross section 
for a^i = 1 while the upper bands (magenta) denotes the cross section for a/^ = 1. The shaded 
regions denote the preferred DM mass region by the CDMS II positive signal. 

We hence conclude that it is very challenging to directly detect the DM cp signal through 
those processes at the LHC. 



V. VECTOR DARK MATTER Z 



The DM candidate can also be a vector boson. In this work we consider a simple extension 
of the SM electroweak gauge group, i.e. adding an Abelian gauge group U{l)x to the SM 
gauge group, extending it to SU{2) l x U{1)y x U (l)^. Note that the SU{2) l x f/(l)y x U{1)x 
is only an effective theory below the scale A where other fields which interact with both the 
SM and the U{l)x in the underlying theory decouple. The Abelian vector field for the 
U{l)x is denoted as Z with a field strength C^^, = df^^Zy — d^Z^. We further assume that 
the quarks, leptons and the Higgs field of the SM do not carry U{l)x quantum numbers, 
and the field Z does not carry quantum numbers of the SM gauge group so that it can be 
the DM candidate. 

We can write down the following dim-6 effective operator 

O^z = \ (0+0) C,,C^\ 
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FIG. 17: Feynman diagrams for ZZ annihilation. Blobs denote the effective vertices induced by 
the dim-6 operators. 

which induces the ZZ — )• hh annihilation via the following vertices 

^(6) _ Oi<l>Z ^ ^f^i, ^(6) _ 0!(f>Z V 

'"ZZhh — ''''''^ii.v^ ) >~'ZZh ~ ' 

In the dimension- six operator, the vector DM Z can either directly annihilate into the Higgs 
boson pair or annihilate via the Higgs-mediated s-channel process into the fermion, Higgs 
boson and vector boson pair; see Fig. [171 It is much simpler as only one operator is present 
now. Needless to say, for a heavy 2, the diagram Fig. ITTT c) is negligible. 



A. Relic abundance 

Since the cross section of ZZ ff is proportional to the fermion mass square, only the 
top quark can contribute significantly to the DM annihilation. When the DM Z is light, say 
80 GeV < m < 150 GeV, it predominantly annihilates into the vector bosons. In the very 
heavy limit, i.e. m ^ rrit, the DM Z will annihilate into WW, ZZ and Higgs boson pairs, 
giving rise to the following simple form of a and b, 

SQa^^m^ 169a^2m2 
" 2567rA4 ' ^ = 20487rA4 ' ^^^^ 

We calculate the relic abundance and plot the allowed parameter space of {mz,A) in Fig. 
M\ for instance, A ~ 700 GeV for mz ~ 100 GeV while A ~ 3000 GeV for mz ~ 2000 GeV. 
Again, since xp barely changes in the whole region, we choose x^? = 25 and obtain an 
interesting relation between m and A as follows: 

0.43 / m \2/TeVV 



cylzi^CDMh^) VlOOGeVV [ A ) ' ^^^^ 
which agrees with the exact scanning results shown in Fig. [TSl But such a simple relation is 
no longer valid when the SM Higgs boson is very heavy because of the large enhancement 
at the threshold uiz ~ Tnh/2. 
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FIG. 18: Allowed parameter set (m^:, A) for a vector DM Z. The upper (lower) boundary of each 
band corresponds to the upper (lower) limit of Q,£)Mh? given in Eq. 1^. 

B. Direct Detection 

Since the vector DM Z only interacts with the SM Higgs boson in our model, it can be 
detected in spin-independent experiments. The relevant Lagrangian is given by 

fe.. = l^^^'^^^"''' ^^^^ 
which gives rise to the effective 2^-nucleon cross sections crf^ and crf^ as follows: 

= 41 « (2.88 X 10- pb) al, {^^) (— ) ) . (46) 

The SI elastic 2^-nucleon scattering cross sections are plotted in Fig. [19] for three choices 
of the SM Higgs boson mass, along with the projected future sensitivity of SuperCDMS 
Stage-A (blue dashed line) and Stage-C (red dashed curve). The cross sections are well 
below all the current bounds from CDMS collaboration, but can be probed at the future 
SuperCDMS experiment at Stage-C. The almost fiat behavior of the black band can be 
easily understood from Eqs. (jS]) and (l46l) as 

Since crfp does not depend upon mz, one can not determine mz from the SuperCDMS 
experiment even if such an excess is indeed observed. But it might be possible to measure 
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FIG. 19: Vector WIMP-nucleon cross sections for spin-independent dark matter search. The shaded 
region is our parameter space which respects the observed WMAP data for a^z = 1: the black 
band for rrih = 120 GeV, the red band for = 300 GeV, and the green band for rrih = 600 GeV. 
See Fig. S] for descriptions of other curves. 

mz via indirect cosmic gamma-ray measurements as discussed below. Needless to say, the 
positive signal of CDMS II would exclude the possibility of vector dark matter. 



C. Indirect search 



The differential cosmic gamma-ray distribution is given by 

'^^ (0.76 X 10~^2g-llp^-2Qgy-l^) J 



^ .5. {^) (I^)%-.-v".. 



We plot the differential distributions as a function of in Fig. I^UT a) for various mz with 
the NFW density profile. For mz ~ 100 — 2000 GeV the distributions are well below the 
current Fermi LAT data. The discovery of the dark matter signal could be possible in the 
high energy region after including an allowed boost factor; see Figs. [20](b) and (c). 

We plot the integrated cosmic photon-fiux as a function of mz in Fig. |2T1 Similar to the 
study of the fermion and scalar DM, we choose two energy thresholds Eth = 1 GeV and Eth = 
50 GeV to mimic the space-based and ground-based telescopes. Assuming J , Afi) /S.VL = 1 
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Log[ E [GeV] ] Log[ E [GeV] ] Log[ E [GeV] ] 

FIG. 20: (a) Predicted gamma-ray spectra for the DM Z with the NFW density profile, where the 
Fermi LAT observation (background) is presented by the red-box (grey-solid) line; (b) Comparison 
between the DM signal plus background and the Fermi LAT observation. Note that a common 
boost factor of 120 is applied to all the DM signals; (c) The comparison between the DM signal 
with the difference of Fermi LAT observation and background. Different boost factors are adopted 
for different DM masses so that the DM signal will not exceed the data. 




fft2 (GeV) 

FIG. 21: Integrated photon flux as a function of mz for energy threshold of 1 GeV (a) and 50 GeV 
(b) (notation of the color bands is the same as Fig. [T8|) . The plot assumes J (^', Ail) Ail = 1; all 
fluxes scale linearly with this parameter. 

and without the boost factor effect, Fermi LAT can probe mz ~ 100 GeV — 1000 GeV while 
MAGIC can probe mz ~ 250 GeV - 1000 GeV. 
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D. Collider search 

The detection of the vector Z can be probed in the VBF process at the LHC as Z only 
interacts with the SM Higgs boson. The production process is given by 

qq — )■ qqVV — > qqZZ. 

The cross section of the production is plotted in Fig. which is too small to be detected 
at the LHC. 

VI. CONCLUSION 

Despite its great success, the Standard Model of particle physics needs to be extended to 
explain the dark matter observation. So far many new physics models have been introduced 
to address DM physics, but our knowledge of the DM candidate is still subject to cosmolog- 
ical observations. For example, the properties of the DM candidate, e.g., its mass and spin, 
remain mysterious. In this work we present a model-independent study of DM physics using 
effective field theory. We add to the SM a new DM field D whose stability is guaranteed 
by a discrete Z2 symmetry. The interactions between the dark matter and the SM fields 
are assumed to be induced by other heavy particles which decouple at the scale A. After 
writing down the subset of dim-6 operators (named dark operators) we are interested in. 
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FIG. 22: Cross section of the vector DM Z pair in the VBF process. 
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we examine the constraints on those dark operators from the rehc abundance measurements 
and direct/indirect searches. The discovery potential for those dark operators at the LHC 
is also studied. 

The dark matter candidate can be a scalar (p^ a fermion % or a vector Z. We found that the 
fermion DM annihilation can be described by three independent operators (C^/^, ^fx^ ^<t>x)i 
the scalar DM annihilation by four independent operators (0</,i, 0</,2, O^f^^, ^f(t>) the 
vector DM annihilation by one operator O^z- Owing to the different setups of various 
experiments, we can probe those dark operators separately. 

Fermion DM x- A positive signal at the CDMS II in the near future will indcate that the 
dark matter he between [1000, 2000] TeV for the operators a)'^ while [900, 1700] GeV for 
the operator a^^. Each individual fermion DM operator could satisfy the cosmic gamma- 
ray spectrum observed by the Fermi LAT since the contribution from the dark matter is 
far below the data. The allowed boost factors can be derived Galactic background agrees 
with the Fermi LAT measurement. For m^^ = 30 GeV ~ 1000 GeV the allowed boost factor 
S = 20 ~ 600. It is possible to detect the heavy dark matter signal from cosmic gamma-ray 
in the higher energy region due to the fact that the dark matter annihilation (including a 
boost factor) produces much more gamma-ray than the background, which is also true in 
the cases of scalar dark matter and vector dark matter discussed below. Furthermore, the 
fermion DM operator might be probed at the LHC via the mono-photon plus missing energy 
and mono- jet plus missing energy signatures. 

Scalar DM (p: The family universal operators are almost ruled out by the CDMS II 
exclusion limit, while the family dependent operators is still allowed with a light scalar 
dark matter. The allowed boost factor is about 25 ~ 300 for the scalar dark matter rri:^ = 
45 GeV ~ 800 GeV), if only af^ is non-zero. For non- vanishing a^i and non- vanishing q;02 
(or q;<^3 = 1, or q;<^i = q;<^2 = Q;</.3 = Q!/</. — 1), the constraints are weaker, i.e. the allowed 
boost factor can be much larger. It is very challenging to search for such scalar dark matter 
at the LHC. 

Vector DM Z : There is no constraint from the CDMS II on vector DM operator as the 
scattering is mediated by the SM Higgs boson which is very small. The allowed boost factor 
is 120 ~ 2000 for the dark matter = 100 GeV ~ 1000 GeV. The cross section of vector 
dark matter signal is too small to be detected at the LHC. 
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TABLE IV: Constraints and potential readies of the dark operators of fermion and vector dark matter in various experiments. Sliorthand 
notations: CII= CDMS II ; Clip = projected CDMS II; SC(A/C)=SuperCDMS(Stage-A/C); B: boost factor need to fit Fermi LAT; {^^): 
full sensitivity; {^/)•. good sensitivity; (— ): not applicable or too small; (x): no constraint; (xx): ruled out. DM mass and scale A are in 
the unit of GeV. 
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TABLE V: Constraints and potential reaches of the dark operators of scalar dark matter in various experiments. Shorthand notations: CII= 
CDMS II ; Clip = projected CDMS II; SC(A/C)=SuperCDMS(Stage-A/C); B: boost factor need to fit Fermi LAT; (W): full sensitivity; 
(^/)■. good sensitivity; (— ): not applicable or too small; (x): no constraint; (xx): ruled out. DM mass and scale A are in the unit of GeV. 
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Appendix A: Relic abundance 



In the computation of the DM (x) rehc density, one assumes that x was in thermal 
equihbrium with the SM particles in the early universe and decouple when it was non- 

1 /2 

relativistic. Once T = < H {H = {Snp/SMpiY'^ is the Hubble expansion rate), 

X stopped annihilating and fell out of equilibrium, and its density remains intact till now. 
The number density is governed by the Boltzmann equation and the law of entropy 
conservation: 



—n^ = -3Hn^ - {(rv)g^ (nj - n^ J - {(rv)^^ {n^n^ - n^,egn^,eq) + Cr, (Al) 



where (av) is the thermally averaged annihilation cross section times the relative velocity, 
the subscripts SA and CA on (av) denote self-annihilation and co-annihilation (annihilation 
with another species 0) respectively, n^^eq denotes the equilibrium number density, and Cr 
is the contribution due to the decay of heavier particles into x- Finally, s is the entropy 
density. For a massive cold dark matter candidate, its equilibrium density is given by the 
non- relativistic limit: 



where m is the mass of the particle species in question. The physics of Eq. (lAll) is as 
following: At early times, when the temperature was higher than the mass of the particle, 
the number density was neq oc T^, x annihilated with its own or other particle species into 
lighter states and vice verse. As the temperature decreased below the mass, dropped 
exponentially as indicated in Eq. ( \A3\i and the annihilation rate T = n (av) dropped below 
H. The X can no longer annihilate and its density per co-moving volume remains fixed. The 
temperature at which the particle decouples from the thermal bath is denoted Tp (freeze-out 
temperature) and roughly corresponds to the time when F is of the same order as H. 

It is usually useful to scale out the effects of the expansion of the universe by considering 
the evolution of the number of particles in a co-moving volume. This is done by using the 
entropy density, s, as a fiducial quantity, and by defining as the dependent variable, Y = n/s 
with Yeg = n^^eq/s- In this case, Eq. ( \Al\i can be rewritten as 



ds 
dt 



3Hs 



(A2) 




(A3) 



dY 

dx 



3Hdx 



1 ds 



3H dx 



1 ds 



(A4) 
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where x = m/T. In radiation domination era, the entropy, as function of the temperature, 
is given by 

which is deduced from the fact that s = {p + P) /T and g^,s is the effective degrees of freedom 
for the entropy density. Therefore, one finds 

? = --■ (A5) 
ax X 

Hence, the Boltzmann equation for the DM number density is obtained 

As is well known, (av) is well approximated by a non-relativistic expansion (obtained by 
replacing the square of the energy in the center of mass frame by s = 4m^ + m^f^): 

{av) =a + b{v^) + 0{{v^)) ^a + Q-. (A7) 

The freeze-out temperature is defined by solving the following equation 

, / , , [45 g mMpi(a + 6b/xF)\ , 

where c is a constant of order one determined by matching the later-time and early-time 
solutions. The result does not depend dramatically on the precise value of c which we will 
choose the usual value c = 1/2. 



Appendix B: Dark matter annihilation 
1. Fermionic dark matter annihilation 

The relevant cross sections for pair of x to annihilate have final states into fermions, 
into Higgs bosons, or into pair of Higgs boson and Z-boson. The Feynman diagrams of the 
annihilations are shown in Fig.O Summing/averaging over final/initial spins and integrating 
over the phase space of the final state particles, we obtain the cross sections of xx annihilation 
shown below. 
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a- XX ff 



The annihilation cross section into a fermion pair is given by 

{xx ff) 

Ncf3f 



127rs/3A4 



{9I + 9I) [s^ -s{m^ + mj) + 4m^m}] + GgigR {s - 2m^) m) \, (Bl) 



where /3 = a/1 — 4m^/ s , /3f = -^1 — Am'j/s and the factor Nc sums over the different color 

combinations allowed in the final state, Nc = 3 for quarks and Nc = 1 for leptons. The 

couplings gi and gn are given in Eq. [71 Expanding avrei in powers of the relative speed 

between the x fermions, Vrei, gives 
2 

m 

Nc <! {gl + gl) (2r) - lOrJ + ll) + 8(7^(7^ rj (3 - 2rJ) j>, (B3) 



« = NcTT-T-Jl - r) [gl + gl + rj g^gR) , (B2) 



487rA4Jl - r) 



where rj = rrif/m. 



b. xx^z^ff 

The annihilation cross section into a fermion pair via the Z-mediated s-channel process 
is given by 

-rUx^Z^ ff) = Nc-^ (gl, + (B4) 

where gzi/R denotes the SM couplings of the Z-boson to fermions and niz (Tz) denotes 
the mass (width) of the Z-boson. We have ignored the masses of the SM fermions in the 
annihilation cross section of the xX ~^ ^ ^ ff process because it only contributes around 
the Z-pole region, i.e. ttidm ~ ^^z, where the masses of the SM lepton and light quarks 
are negligible. The corresponding leading terms in the non-relativistic expansion are 

SttA* (4m2 - m|) + m|r| 



i = -Nc<^, (Jk + le^' + ^O'"!'/ (B6) 
^mTiA'^"^'' "^"l (4m2 - m|)' + m|r| ^ 

As to be shown later, this annihilation channel is dominant around the Z-pole region, i.e. 
m ~ \mz^ but its contribution decreases so rapidly, once goes beyond the Z-pole, that 
it is overwhelmed by other annihilation channels. 
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c- XX hZ 

The annihilation cross section into a pair of Z-boson and Higgs scalar shown in Fig. Wis) 
is given by 

o- (XX hZ) 

,2 



A (s, m|, ml) {{s + 2m^) A (s, m|, m^) + 12s m| (3m^ - -s) ^ (B7) 



967rs3/3A4 
where the function A being defined as 

A (x, y,z) = {x - y - zf - Ayz. 

The non-relativistic expansion gives rise to the following leading terms 

2 2 
2 2 

-64 (3r^ - 16rM - 16r|) + 8 (iSr^ - 2r|) (r^ - r|) - 11 (r^ - r|)' . (B9) 
where r/,(z) = mh{z)/m. 

d- XX ^ Z ^ hZ 

The annihilation into a pair of Z-boson and Higgs scalar can also be induced by the 
Z-mediated s-channel process, see Fig. [2](d). This annihilation channel opens only if 277i > 
777^ + TUh, which implies s ~ 4777^ > 4rr7|, thus its contribution is highly suppressed so that 
we can ignore it in this study. 
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e- XX 



w+w- 



For this channel, 



levrA^ 



17 



4 

145 
"28" 



1 - 



-gM9 + 55r^-^rt. + 6r^ 



1-^ 



(BIO) 



(Bll) 



where rw{z) = mw{z)lm. 



2. Scalar dark matter annihilation 

When the scalar is the DM, we must consider a large number of annihilation processes, 
e.g. ff/WW/ZZ/hh and ipip ^ h ^ ff/WW/ZZ/hh, see Fig. M The latter, 

Higgs-mediated s-channel processes, contributes significantly around the Higgs resonance 
region only. 



a. ipip ff 



In the scenario-A, the annihilation cross sections of ipip — )■ // read as 



- - ff) - 'Vc^/^f + NcM.f" ^ 



167r/3A 



327r/3A 



(B12) 



{s - mlf + mlTl' 

where the first and second term corresponds to Fig. [TOT a) and (b), respectively. Expanding 
avrei in powers of the relative speed Vrei gives the following leading terms 



V 



c 



167rA4 

2 



(1 ' '■?) 



2\3/2 



2 2 

mjv 



N, 



V 



c 



2 



1-rj 



{Am? -ml) +mlVl 
2a), (2rJ + l) 



(B13) 



I 2 2 



[-48m6 + 8 (I2m^ + ml) m^ + ml [ml - ?>2m)) m^ + 2mjml] Lg^^^ 
{{Am^ - mlf + mlTl)' J 
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Similarly, in the scenario-B, the leading terms are given by 



V 



c 



2 



(1 - -'?) 



2\3/2 



2a 



f 

f<t> I ^ 



+ a 



2 2 

nijV 



01 



{Am? - mlf + mlTl 



(B15) 



V 



c 



1287rA4 

2 



V 



I 2 2 



6. (^v9 ^ ly VF/zz 

The annihilation cross sections of (pip annihilation into vector bosons are given by 

a [(pip -> VV) = 6v — —. + -f^- — (B17) 



647rs/3A^ 



s-ml) +mlTl 



where the subscript V denotes the type of vector boson and f3y = a/1 — Amy/s. 6v counts 
for the identical particle: 6w = 1 and 6z = 1/2. The first and second term in Eq. IB 171 
corresponds to Fig. [TOlfc) and (d), respectively. The non-relativistic expansion gives rise to 
the following leading terms 

5v 



ay 



327rm3A4 



m^ — mt 



X < a 



</)3 



(4m^ — Amlm^ + Sm^) + 



2 4 



4m'* — Am'^m'^ + 3my 



[m 



I - Am^f + min 



(B18) 



X 



2567rm3A4 ^/^^^ niy 
f (4m^ — m^) 



(^(m2 - 4m2)^ + m^r^ 



16m^ + 4 (3m^ - 20m^) - 4 (3m^m^ - 31mt^) 



/I/O o \ o 

+3my (^m^ + 2Amyj m + 6m 



(B19) 
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c. (^9? — )• hh 



The cross section of (p(p hh reads as 



where 



The first and second term in Eq. IB20I corresponds to Fig. [TOlfe) and (f), respectively. The 
non-relativistic expansion gives the following leading terms 



m? — m\ 



n 2 4 4 

9 V 



^ ' {Am^ - mlf + mlT^ 



^ (4c2m^ + Cit;^) [Ac2m^ {Sm"^ — 2m\) — Cif ^ (m^ — 2m^)] 



(B21) 



5127rm^A'^ A^m^ — 
^ 9a^im^ (-SOm^ + 120m^m^ - 33mlm^ + 2ml) 



(B22) 



3. Vector dark matter annihilation 

The vector DM Z can either directly annihilate into the Higgs boson pair or annihilate 
into the fermion, Higgs scalar and vector boson pair through the Higgs-mediated s-channel 
process, see Fig. [17] 

a. ZZ — )• hh 

The annihilation cross section of ZZ — )■ hh is given by 

which gives rise to the leading term in the non-relativistic expansion as follows: 

3«L^^ / (13m2 - 10m2) ^ ^ 

2567rA4 V 20487rA4 71^72 ^ ^ 
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b. ZZ ^ h ^ hh 



The annihilation cross section of ZZ pair into the Higgs scalar pair via the Higgs- mediated 
s-channel process is given by 

which yields the following leading non-relativistic expansion terms 



5127rA4 



c. ZZ^h^ ff 



The annihilation cross section of h ^ ff is given by 

(7 (ZZ ^ /i ^ //) = ATg ^ ^ i, (B28) 

which yields 

a = Nc !\J ^ , (B29) 

We note that both a and 6 are proportional to m/ so that they are negligible in the limit of 
m S> mf. 



d. ZZ^h^ WWjZZ 

The annihilation cross section of ZZ into vector boson pair is given by 

^(77^h^VV\ S "^^^^^ {6m^-Asm'' + s^){12ml-Asml. + s^) 

a{ZZ ^ VV) ^dv-—TTj ^ — , (B31j 
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where 5w = ^ and Sz = 1/2. The non-relativistic expansion gives rise to the two leading 
terms as follows: 



ay = Sv \ \-rl— ' ^ (B32) 



5^ "' ^ ""^^ 208m8 - 4 (25m2 + eSm^.) 



+4 (I9mt. + Almlml) + (24m^ - 103m^m^) + 30m^m^ [> (B33) 



Appendix C: Dark matter direct detection 



Fermionic DM-nucleon interaction 



When the DM is a fermion x, its coupling to quark, leading to elastics scattering from a 
nucleus, can be parametrized as Eq. [T31 of which the first term results in a spin-independent 
scattering from a nucleus while the second term leads to a spin-dependent scattering from 
a nucleus. 

Consider the SI scattering first. When a WIMP interact with quarks via a vector- like 
interaction, e.g. given by 

f^v = bqXl/^Xqi^Q, (CI) 

where bg is a shorthand notation of the WIMP-quark vector coupling given in the first term 
of Eq. [131 In this case, the contributions of each quark in the nucleus add coherently and 
large cross sections result for large nuclei. The WIMP-nucleus cross section in this case is 

- ~( \ 72' y^^) 

where = 2Z bp + {A — Z) bn with bp = 26„ + bd and bn = b^ + 2bd- It is convenient 
to consider the cross section with the single nucleon for comparing with the experiment. 
Through a simple algebra calculation, we obtained the cross sections of WIMP-proton and 
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WIMP-neutron as given below: 



SI _ ^"dm^I 1 
'''''' ~ 7r(m,,M + mp)'l6A4 



1 2 



(6.98 X 10-5 pb) ( - + 0.013«^^) , (C3) 



2 2 



«-<-0.162a^^)'. (C4) 

Now consider the SD scattering. WIMPs could also couple to the spin of the target 
nucleus through an "axial-vector" interaction 

>Ca = dqXl^L75XQl^l5q, (C5) 

where dg is the shorthand notation of the axial coupling given in the second term of Eq. 
The spin-dependent WIMP-nuclei elastic scattering cross section can be expressed as 



a. 



SD ^ ^^iii'DM"''N 



where J is the total angular momentum of the nucleus. Aat (oc 1/J) depends on the axial 
couplings of WIMPs to the quarks, 

_ ap {Sp) + an {Sn) 
An = J . 

where 

ap = d^Al + d,A2 + 4A^, a„ = d^Al + rf^A^ + ^A^. 

Here, A's are the fraction of the nucleon spin carried by a given quark. Their values are 
measured to be [tsI 

AP = A;J = 0.78 ± 0.02, A^ = A;^ = -0.48 ± 0.02, A^ = A^ = -0.15 ± 0.02. 

{Sp) and {Sn) are the expectation values of the total spin of protons and neutrons, respec- 
tively, whose values depend on the number of protons and neutrons in the nucleus being 
considered, namely, whether it is odd or even. For odd-proton nuclei the spin-dependent 
WIMP-nucleus cross section is mainly due to the WIMP-proton interactions, whereas for 
odd-neutron nuclei it is dominated by WIMP-neutron scattering. For even-even nuclei the 
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spin-dependent cross-section is highly suppressed. For a proton and a neutron as a tar- 
get, Eq. [T7]is transformed into the cross section form WIMP-proton (neutron) interactions 
with the proton (neutron) spins {Sp^n) = 1/2 and J = 1/2. Thus, the SD cross section of 
WIMP-proton and WIMP-neutron are given by 

(rp -t r\ 4 
(«J^ + «J'^-4.53a^;,)', 

2. Scalar DM-nucleon interaction 

When the DM is a scalar ip, it can be detected in the spin-independent experiments via 
the scalar interaction with nucleus 

which leads to the scattering amplitude 

{M) = C^, {qq) , 

where ( ) denotes an average and sum over the spins of the initial and final state quarks, 
respectively, and the coefficient C^q is given by 

_ 1 f af4,v , g^it;^ nig 

The matrix element {qq) of quarks in a nucleon state is given in jl6| by 

{qq) = ^^/^^^'"^ (light quarks) ; {qq) = ^^^^^Jt^^ (heavy quarks) . 
Summing over quark flavors, we obtain the <y9-nucleon couplings: 

q=u,d,s ^ q=c,b,t ^ 
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where [79|, |S 

c\ non —I— n nn/i n c\^f\ —I— Ci nn 



/}^^ ^ 0.020 ± 0.004, /^^^ ^ 0.026 ± 0.005, /^^ ^ 0.118 ± 0.062, 

^ 0.014 ± 0.003, ^ 0.036 ± 0.008, /^"^ ^ 0.118 ± 0.062 . (C9) 



The first term in Eq. IC8I corresponds to interaction with the quarks in the target nuclei, 
while the second term corresponds to interactions with the gluons in the target through a 
quark loop diagram. f!^J is given by 1 — f!^J — f!j?J — f!^J ^ 0.84, and analogously, /jf"* ~ 0.83. 
Note that we have to sum over the couplings to each nucleon before squaring because the 
wavelength associated with the momentum transfer is comparable to or larger than the 
size of nucleus, the so-called "coherence effect". Due to such a coherence effect with the 
entire nucleus, the cross section for spin-independent interactions scales approximately as 



8l| 



the square of the atomic mass of the target nucleus 

The total spin-independent WIMP-nuclei cross section at zero momentum transfer is 
given by 

4n = , , [Zf^'J + {A-Z) /(^)] ' , (CIO) 

where is the target nuclei's mass, and Z and A are the atomic number and atomic mass 
of the nucleus. In order to compare with the experimental sensitivities and limits which are 
often described in terms of the dark matter elastic scattering with nucleons, we derive the 
V9-nucleon cross section as following: 

2 2 

SI _ r.(p)]2 SI _ [fWr 



4vr {mnM + m,f ^^^^ ^ ' An {m^M + m^f ^^^^ ^ " 

In practices, we found a^p ~ o"^^. 

3. Vector DM-nucleon interaction 

Since the vector DM Z only talks to the Higgs boson in our model , it can be detected 
in the spin-independent experiments. The relevant Lagrangian is given by 

J-ZZqg - 

which gives rise to the following scattering amplitude 

^Mzz,, = [{Pi ■ Ps) g'" - KK] el (P3) e, {p,) q {p,) q {p,) . 

A 
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In the extreme non-relativistic limit, Pi = P3 = {mr)M,0) and the polarization vector of 
the heavy vector Z are purely spatial, £^ (^1,3) = (0,61^3). The above amplitude is thus 
simplified as 

iMzZqq = Czq (ps) (Pl) Q (Pa) Q fe) , 

where 



mjj nig 
Is? m\ 



The effective WIMP-nucleon scattering cross section is thus given by 



SI 

Zp,Zn 



m. 



p,n 



Zitf + {A-Z) /g 



(p,n) 



An {mo + nip^n) 
where the Z-nucleon couplings, fzq"\ are given by 

r{p,n) _ Sr~^ Ap,'n-)p f^p,n 2 „(p,n,) ^ nip^n 

JZq — ?, JTa ^Zq 1" T^JTg 2. " 

q=u,d,s ^ q=c,b,t ^ 

Again, using rup ~ m„ and m£, ^ fTip^n, we obtain the effective WIMP-nucleon cross section, 
^zpi= (^zi) as following: 



SI 
Zp 



a 



4>z 



m 



D 



mi 



1 2 



(2.88 X 10"^° pb) 



rriD 



TeV 



100 GeV 



V 100 GeV/ \ A J \ nih 

The SI elastic 2^-nucleon scattering cross sections are plotted in Fig. ^\ along with 
the projected future sensitivity of SuperCDMS, Stage-A (blue dashed line) and Stage-C 
(black dashed curve). The cross sections are well below all the current bound from CDMS 
collaboration, the can be probed at the future experiment by SuperCDMS at Stage-C. 



Appendix D: Dark matter indirect detection from gamma-ray 

WIMPs annihilation in the halo may lead to a flux of gamma-rays, with both continuum 
and line contributions. The former can be produced as continuum photons from flnal state 
radiation and the cascades of other annihilation products, whereas the latter can be emission 
from loop-diagrams to 77, 7Z or 7/1 flnal states. Observation of monochromatic gamma rays 
would provide a "smoking-gun" signal for the existence of WIMPs in the halo. 
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The produced gamma ray flux from the annihilation of dark matter particles can be 
expressed as 

= / p^di (Dl) 

dVt dE ^ dE^^ ' AttttiIm Ji.o.s 
where p is the dark matter density proflle, (o"jf ) and dN^/dE^ are, respectively, the thermally 
averaged annihilation cross section times the relative velocity v and the differential gamma 
spectrum per annihilation coming from the decay of annihilation products of flnal state i. 
The integration is taken along the line of sight. 

All the halo model dependence isolated in the integral which is given in the dimensionless 
form Q 

J(^,An) = — \ — ( — -) [ dn [ p'^dl, 

^ ' ^ 8.5kpc V0.3GeV/cm3; An J^^ Ao./ 

where \l/ is the angle away from the direction of the galactic center that is observed while 



AQ the solid angel of the fleld of view centered on \1/ = 0. A 



a spherical solid angle AQ, one obtains the photon flux as 82 1 



ter averaging J (\E', AQ) over 



dE^ 



where the function J contains all information about the dark matter distribution in the 
halo. The above equation can be factorized into two parts: J (\E', AQ) AQ (astrophysics) and 
other terms describing the dark matter annihilation and the fragmentation of its annihilation 
products (particle physics). The latter can be determined from microscopic measurement, 
but the former exhibits considerable uncertainties involved in the distribution of dark matter 
in the Galactic center region. For example, at AQ = 10~^ sr, characteristic of ground- 



based Atmosp^ 
NFW proflle 



leric Cerenkov Telescopes (ACTs), typical values of J range from 10'^ for the 
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66| . to about 10^ for the proflle of Moore, et at jSSj. Furthermore, as 



the distribution of dark matter might be clumpy, this inhomogeneity effect would enhance 
indirect detection rates by a "boost factor" , deflned as 

B^ri- (D3) 

If the dark matter were locally distributed completely evenly, the boost factor would be 
equal to one. Small-scale structure of dark matter, however, enhance this quantity to a 
large value. 
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The particle physics model dependence enters through all the other factor in Eq. ID1[ 
The energy integral is roughly J dE dN^/dE ~ 0.5 for all i, but the energy distribution 
depends significantly on the annihilation channel. Note that the resulting photon spectra 
depend only on the initial energies of the annihilation products and not on the details of 
the WIMPs annihilation process. The spectra have been studied using PYTHIA 8J| , and a 



simple analytic fit has been presented in Refs. 82|, |85|] for the most important annihilation 
channels as follows: 

dN^ ^ _a_^-bx 
dx x^-^ ' 

where x = E^/moM and (a, 6) = (0.73,7.76) for WW and ZZ, (1.0, 10.7) for hi, (1.1, 15.1) 
for tt, and (0.95, 6.5) for uu. In most cases the spectrum produced does not vary much. The 
only exception to this is the somewhat harder spectrum generated through annihilation to 
r+r^. Particle physics models, however, very rarely predict a dark matter candidate that 
annihilates mostly to r+r^. In this work we focus our attention only on the annihilation 
channel involving the heavy fermion, gauge boson and Higgs boson. 



85| 



The integrated photon flux above some photon energy threshold Et^ is 

%{Eth) = (5.5 X 10-1 V^cm-^) (^M^^ J (v];, Afi) Afi 



X 



where the sum is over all possible annihilation channels. Detectors also have upper cutoffs, 
but these are typically irrelevant, as the energy distribution falls steeply with energy. In 
this work we consider two representative Eth'- 1 GeV, accessible to space-based detectors, 
and 50 GeV, characteristic of ground-based telescopes. 
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